Astronomy & Astrophysics manuscript no. banl2059 © ESO 2009 

May 7, 2009 



Letter to the Editor 



Propagating waves in polar coronal holes as seen by SUMER & EIS 

D. Banerjee 1 , L. Teriaca 2 , G. R. Gupta 1 ' 3 , S. Imada 4 , G Stenborg 5 , and S. K. Solanki 2 * 6 



On 
O 
O 
<N 

>> 



O 



> 



Indian Institute of Astrophysics, Bangalore 560034, India 

Max-Planck-Institut fur Sonnensystemforschung (MPS), 37191 Katlenburg-Lindau, Germany 
Joint Astronomy Programme, Indian Institute of Science, Bangalore 560012, India 
National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 
Interferometrics, Inc., Herndon, VA 20171, USA 

School of Space Research, Kyung Hee University, Yongin, Gyeonggi 446-701, Korea 



ABSTRACT 



Context. To study the dynamics of coronal holes and the role of waves in the acceleration of the solar wind, spectral observations 
were performed over polar coronal hole regions with the SUMER spectrometer on SoHO and the EIS spectrometer on Hinode. 
Aims. Using these observations, we aim to detect the presence of propagating waves in the corona and to study their properties. 
Methods. The observations analysed here consist of SUMER spectra of the Ne vm 770 A line (T = 0.6 MK) and EIS slot images in 
the Fe xn 195 A line (T = 1.3 MK). Using the wavelet technique, we study line radiance oscillations at different heights from the limb 
in the polar coronal hole regions. 

Results. We detect the presence of long period oscillations with periods of 10 to 30 min in polar coronal holes. The oscillations 
have an amplitude of a few percent in radiance and are not detectable in line-of- sight velocity. From the time distance maps we find 
evidence for propagating velocities from 75 km s _1 (Ne vm) to 125 km s _1 (Fe xn). These velocities are subsonic and roughly in the 
same ratio as the respective sound speeds. 

Conclusions. We interpret the observed propagating oscillations in terms of slow magneto- acoustic waves. These waves can be 
important for the acceleration of the fast solar wind. 
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1. Introduction 2. Observations 



* Propagating radiance oscilla tions were detected in polar plumes, The data analysed here were obtained on 8* h and 15* April 2007 

first by lOfman et al.l {l 997) using UVCS/ SoHQ and later by during a Hinode/SUMER joint observing campaign as part of 

^ I [ Deforest & Gurmaddl998h withEIT/SoHQ. lQftrian et al.l(ll999L the Hinode Observing Programme (HOP) 45/Joint Observing 

2000a) identified the observed radiance oscillations as propagat- program (JOP) 196. They consist of time series taken in the 

• • \ ing slow magneto-acoustic waves. A number of studies using the south polar coronal hole by the Solar Ultraviolet Me asurements 

^ ■ CDS/SoHO spectrometer have reported oscillat ions in the po- Of Emitted Radiation (SUMER. IWilhelm et al.ll 19951) spectrom- 

' lar coronal hol es up to 25 " a bove the lim b (e.g.. iBanerjee et al.l eter aboard the Solar and Heliospheric Obs ervatory (SoHO) and 

2QQQ. l2001allb |; IO , Shea et al.ll2006ll2007h . lPopescu et al.1 (12005b . by the EUV i maging spectromete r (EIS, ICulhane et aD l2007h 

using SUMER/SoHO, detected radiance fluctuations with peri- aboard Hinode dKosugi et al.ll2007h . For SUMER, the 1" x 120" 

ods from 10 to 90 min up to 15" above the limb. These studies slit was centred on the limb and spectra were acquired from 

point to the presence of compressional waves, thought to be slow 19:13 to 20:47 UTC on % th and from 10:44 to 14:45 UTC on 

magneto-acoustic waves. In this letter, for the first time to our 15 r \ with a cadence of 18.12 s. For EIS, the 40" wide slot was 

knowledge, simultaneous use of the SUMER and EIS/Hinode used to obtain 40" x 160" images in several spectral lines over 

spectrometer were used to study these propagating disturbances the time interval 18:42 to 20:58 UTC on % th and 10:54 to 15:57 

in the off-limb regions of the polar coronal holes. We construct UTC on 15^ April. The EIS data consist of a series of elemen- 

time distance maps to study the properties of wave propagation tary rasters each formed by two slot images displaced by 20" in 

and use wavelet analysis to establish their periods. Spectroscopic the X direction to maximise the chances of overlapping with the 

observations have the advantage of a narrow temperature re- other instruments. Each slot image was exposed for about 7.5 s. 

sponse (by isolating specific spectral lines) and of allowing the As a result, for each dataset, we have two time series with a ca- 

study of resolved and unresolved plasma motions by measuring dence of 19.3 s that we identify hereafter with sloto (East) and 

the Doppler shift and width of the observed profiles. These ob- sloti (West). Before the start of the temporal series, raster images 

servables provide important constraints in establishing the nature were obtained with SUMER and EIS to allow the co-alignment 

of the observed oscillations. of the different instruments. 
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Fig. 1. Left (right) panel shows the location of different slits and 
slots on the Fe xn 195 A EIS context raster on 8 th April (15 th 
April). The context rasters are made of five adjacent slot im- 
ages, covering a 200" x 320" area. In both panels, radiance 
iso-contours (white) are from EIT Fe xn 195 A images. The two 
white rectangular boxes are the location of the EIS slots (solid 
sloto, dashed sloti). The solid vertical blue line represents the lo- 
cation of the SUMER 1" wide slit whereas the yellow dash-dot 
elongated box gives the location of the 4" wide CDS slit. 



All the data were reduced and calibrated with the standard 
procedures given in the SolarSoft (SSWfl library. SUMER data 
were decompressed, corrected for response inhomogeneities 
(flat-field) and for geometrical distortion (de-stretch). The data 
series of 8 th April was analysed only until 20:47 UTC due to 
spurious flaring at the edges of the detector occurring after this 
time. For EIS, after applying the standard reduction and calibra- 
tion provided by the eis_prep procedure, data were corrected for 
the spacecraft jitter in the Y-direction (the jitter in the X direc- 
tion is less than 1" and can be neglected) by using housekeep- 
ing data. Finally, the movement of the slot image on the detector 
due to thermal variations during the orbit was corrected. The dis- 
placement in the dispersion (Solar X) direction was obtained by 
measuring the position of the edge of the Fe xn 195 slot image 
over time. The displacement in the Y direction is taken equal 
to 1.5 times that in the X direction (Imada, 2009 in preparation). 
The validity of the latter assumption was verified by checking the 
limb Y position vs time. All images from Hinode were converted 
to SoHO view (LI) and co-aligned in two steps. First, the in- 
ternal offset between the long-wavelength and short- wavelength 
EIS CCDs and the wavelength dependent inclination of the spec- 
trum were obtained (Kamio & Hara, private communication) 
and accounted for. Then, choosing SUMER as the reference, the 
EIS rasters were aligned with the SUMER rasters. We estimate 
the alignment to be accurate within 2". Fig. [T] shows the loca- 
tion of the slits and slots of the different instruments. Here the 
white rectangular boxes represent EIS sloto and sloti . The ver- 
tical lines indicates the location of the SUMER (solid blue) and 
CDS (dash-dotted yellow) slits. Due to the marginal overlap with 
the other instruments and the low signal in the off-limb spectra, 

1 http : //sohowww . nascom . nasa . gov/solar sof t/ 
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Fig. 2. Enhanced (see text) maps of radiance variation along the 
slit (Solar Y direction) with time for Ne vm as recorded by 
SUMER (top panel) and Fe xn 195 A as recorded by EIS (bot- 
tom panel) on 8 th April 2007. The vertical black dashed line on 
the SUMER enhanced radiance map depicts the starting point 
of the EIS time series (shown in the bottom panel). The slanted 
dotted yellow lines correspond to disturbances propagating with 
a speed of 75 km s" 1 and a period of about 15 minutes, as deter- 
mined from the SUMER data. 



3. Results 

Maps of the radiance along the slit vs time (xt slices) were first 
built using the SUMER Ne vm integrated line radiances and EIS 
Fe xii radiances averaged over 5" in the X direction at the posi- 
tion overlapping the SUMER slit. The resulting maps were then 
contrast enhanced and low-pass filtered to produce the maps 
shown in Figs. [2] and [3] for the 8 th and 15 th April datasets, re- 
spectively. EIS radiance maps were normalised by dividing by 
the average along the Y direction (to remove the strong limb 
signature), the result elevated to the power of 3 (to increase the 
contrast between high/low radiance regions), and, finally, low- 
pass filtered. SUMER enhanced radiance maps were obtained 
as the median over 3x3 pixels of the quantity exp(a + 0.15 x b) 2 
where b is the original image, and a is a highly smoothed version 
of the original image obtained by convolving with an isotropic 
kernel. The difference in the treatment of the data from the two 
instruments depends on the different noise level and image scale. 
In Fig. O we note that the EIS time series starts at a later time 
(denoted by a dashed vertical line in the top panel of the figure). 
On xt slices, the presence of alternate bright and dark regions 
indicates the presence of oscillations. Moreover, diagonal radi- 
ance enhancements (marked by dotted yellow lines and dashed 
blue lines in Figs. [2] and [3]) are a signature of propagating distur- 
bances. Thus, from such maps, it is possible to estimate periods 
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Global Period at max. 
power (< 34.0 min.) 
= 14.12 min. 
Prob. level: 99-100% 




! II I I 



Hi- 



Time (minutes) 
FeXII Wavelet 



Global Period at max. 
power (< 20.1 min.) 
= 13.89 min. 
Prob. level: 99-iC0% 



Global Wavelet 




Time (minutes) 

Level of Probability 



20 30 40 50 

Time (minutes) 

Level of Probability 



100 1000 10000 



60 80 



Time (minutes) 



Fig. 4. Wavelet analysis for the 8^ April data at Solar Y~ -995" for both Ne vm 770 A (left) and Fe xn 195 A (right). Both light 
curves are obtained by averaging over 9" along the slit. In each set we show the relative (background trend-removed) radiance (top 
panels), the colour inverted wavelet power spectrum (central panels), the variation of the probability estimate associated with the 
maximum power in the wavelet power spectrum (marked with white lines) (bottom panels), and the global (averaged over time) 
wavelet power spectrum (right middle panels). Above the global wavelet, the period (measured from the maximum power of the 
global wavelet), and the probability estimate, are given. 
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Fig. 5. Wavelet analysis for the 15 th April data at Solar Y~ -995" for both Ne vm 770 A (left) and Fe xn 195 A (right). Both light 
curves are obtained by averaging over 9" along the slit. Panels as in Fig.|4] 



and projected propagation speeds. From the SUMER data on 
both dates we measure a propagation speed of -75 km s" 1 and 
a periodicity of about 15 minutes. The propagation fronts are 
much more visible on the 8 r/ * April data (see top panel of Fig. [2]) 
but are also visible on the 15 th April data. In the case of the EIS 
data, no definitive conclusions can be drawn from the 8 th April 
data (see bottom panel of Fig. O, while propagating fronts are 
clearly visible in the 15 th April data (see bottom panel of Fig. [3]). 
However, both the propagation speeds and the periodicity are 
higher than for the SUMER data, with values of -125 km s" 1 
and 30 minutes, respectively (see the dashed blue lines in the 
bottom panel of Fig. [3]). The observed propagation speeds are ac- 
tually a lower limit as the field lines may form an angle with the 



plane-of-sky. However, these angles are very likely small and the 
speeds are most probably subsonic. For the propagation speeds, 
we estimate the uncertainty to be around 10%. The ratio of the 
Fe xii to Ne vm propagation speeds, about 1.7, is close to the ra- 
tio of the sound speeds of about 1 .5 (at the two formation temper- 
atures of 0.6 and 1.3 MK). The observed radiance perturbations 
have an amplitude of a few percent. No evidence of propagat- 
ing disturbances can be found in the LOS velocity xt maps (not 
shown). 

To better study the properties of the propagating disturbances 
seen in the enhanced-radiance xt slices in Figs. [2] and [3l we 
make use of wavelet analysis and focus on individual locations 
(heights) in the off-limb corona. For this purpose we extracted 
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Fig. 3. Enhanced radiance maps as in Fig. O but for the 15 th 
April data. Black vertical stripes correspond to data gaps in the 
EIS time series. The horizontal white dashed line on the EIS 
enhanced radiance map depicts the upper edge of the SUMER 
FOV. The slanted dotted yellow lines correspond to disturbances 
propagating with a speed of 75 km s" 1 and a period of about 
15 minutes as determined from the SUMER data. The slanted 
dashed blue lines on the EIS map correspond to disturbances 
propagating with a speed of -125 km s" 1 and a period of about 
30 minutes as determined from the EIS data. 



radiance vs time curves at a given height by averaging over 
9"in the Y direction for both SUMER and EIS original (with- 
out filtering or contrast enhancement) radiance maps. In Fig. 01 
(8 th April) and Fig. 0(15^ April), we show representative ex- 
amples of the type of oscillation measured at Y= - 995 " 
(~ 28" above the limb). The top panels of Figs. 01 and \5\ show 
the variation of the radiance (hereafter we will use the term radi- 
ance for trend- subtracted line radiance). Details on the wavelet 
analysis, which provides infor mation on the temporal var iation 
of a signal, are described in iTorrence & Compol (1 19981) . For 
the convolution with the time series in the wavelet transform, 
we ch ose the Morlet function, as defined in ITorrence & Co mpo 
(1998). The light curves shown in the upper panels had their 
background trend removed by subtracting from the original time 
series a 60-point (« 18.1 min) running average for 8 th April 
data and a 100-point (« 30.2 min) for SUMER and 150-point 
(« 45.3 min) for the EIS running average for 15 th April data. 
In the wavelet spectrum (middle-left panels), the cross-hatched 
regions are locations where estimates of oscillation periods be- 
co me unreliable. This is the so-called cone-of-influence (COI, 
see ITorrence & Compo| [T998). As a result of the COI, the maxi- 
mum measurable period is shown by a dotted line in the global 
spectrum plots (middle-right panels). Above the global wavelet 
spectrum of Figs. 0] and [5] we show the prevalent period, mea- 
sured at the location of the maximum of the global wavelet spec- 
trum, together with an estimate of the probability that this oscil- 
lation is not due to noise. The probability estimate was calcu- 
lated using the rand omisation me t hod w ith 200 permutations as 
outlined in detail in lO'Shea et aD d2001h . 



Below the wavelet power spectrum, in the lower panels, we 
show the variation of the probability estimate, calculated us- 
ing the randomisation technique, associated with the maximum 
power at each time in the wavelet power spectrum. The location 
of the maximum power is indicated by the over-plotted white 
lines in the middle-left panels. Figure 0] shows that the oscilla- 
tions detected in the 8 th April data have periodicities with max- 
imum power at ~14.1 min in Ne vm (SUMER) and ~13.9 min 
in Fe xn (EIS). Note that the wavelet technique reveals highly 
significant power in Fe xn (also at other heights), although no 
clear oscillation could be discerned by eye in the EIS enhanced 
radiance map (see bottom panel of Fig. [2]). On the other hand, for 
15 th April data, which have a longer time series, Ne vm shows 
its maximum power around 1 8 min and Fe xn shows maximum 
power around 30 min, consistent with the results from the anal- 
ysis of the xt slices. Also note that there are multiple peaks seen 
in the global wavelet spectra (right panels), which implies that 
the signal could be composed of multiple modes of oscillation. 
To study the oscillation behaviour at different heights we also 
look at the wavelet results at Y= -1003" and plot the results 
in Figs. [6l and [71 respectively (available online only). For the 8 th 
April data at Y= -1003", SUMER reveals a strong peak at 14.2 
min, while EIS shows a peak at 15.15 min. For the 15 th April 
data at Y = -1003" SUMER shows a strong peak at 20 min, 
whereas EIS shows a peak at 30 min. Finally, we also attempted 
to search for oscillations in the LOS velocity (by fitting the pro- 
files obtained using the same binning as for the radiance wavelet 
analysis). No evidence of oscillations with an amplitude larger 
than 2 km s" 1 can be found in the LOS velocity data. This is the 
accuracy achievable with SUMER when observing the Ne vm 
line in the first order of diffraction (at the Signal to Noise values 
of relevance here). 

4. Conclusions 

The first observational detection of longitudinal waves came 
from analysing polarised brightness (density) fluctuations in 
white light data. Fluctuations with periods of about 9 min 
wer e detected in coron al holes at a height of about 1.9Rq 
by lOfman et~ai] (1 19971) using th e white light ch annel of 
UVCS/SoHO. In a follow-up study. lOfinan et alJ (l2000bl) deter- 
mined the fluctuation periods to be in the range of 7 to 10 min. 
The propagation speeds of the fluctuations indicated values in 
the range of 160 to 260 km s" 1 at 2 R@, whic h is slightly slower 
than the acoustic speed at those heights. iDeforest & Gurmanl 
(1998), using EIT 171 A reported detection of outwardly prop- 
agating radiance perturbations at distances of 1.01 to 1.2 Rq, 
gathered in quasi-periodic groups of 3 to 10 periods, with pe- 
riods of about 10 to 15 min. The projected speeds are about 
75 to 150 km s" 1 and the relative amplitude (in density) was 
about 2 to 4 %. Usually these waves are observed propagating 
along the assumed coronal magnetic structures and, thus, along 
the magnetic field. Their speeds are usually much slower than 
the expected coronal Alfven speed, which leads to their inter- 
pretation as longitudinally propagating slow magneto-acoustic 
waves. Slow magneto-acoustic waves follow magnetic field lines 
and propagate at the local sound speed. We detect the presence 
of long period oscillations with periods of 10 to 30 min in po- 
lar coronal holes within the range of 1 to 1.2 Rq, with a clear 
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signature of propagation with velocities from 75 to 125 km s , 
depending on the temperature of line formation. The measured 
propagation speeds are subsonic, indicating that they are slow 
magneto-acoustic in nature, which is consistent with earlier re- 
ports. Note that this detection has been confirmed through the 
analysis of data from two separate spectrometers on-board two 
different satellites. Thus we feel that this simultaneous detec- 
tion makes the result very robust. We also find that the prop- 
agation speed in Fe xn (-125 km s" 1 ) is higher than that in 
Ne vm (~75 km s" 1 ), as shown in the bottom panel of Fig. [3] 
This may be a temperature effect, as the ratio of the Fe xn to 
Ne vm propagation speeds, about 1.7, is close to the ratio of the 
sound speeds. Different propagating speeds observed in differ- 
ent lines may also be interpreted as an indication of the presence 
of structures with different temperatures along the line of sight 
(e.g., weak plumes). The observed region can be either a bundle 
of magnetic threads of different temperatures, or have a trans- 
verse temperature profile. However, with the available data this 
statement remains a conjecture only. Finally, we note that the ob- 
served waves have no detectable signature in the LOS velocity, 
enforcing the idea of compressive longitudinal magneto-acoustic 
waves. 
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Fig. 6. Wavelet analysis for the 8 th April data at Solar Y~ -1003" for both Ne vm 770 A (left) and Fe xn 195 A (right). Both light 
curves are obtained by averaging over 9" along the slit. Panels as in Fig. [H 




Fig. 7. Wavelet analysis for the 15 th April data at Solar Y~ -1003" for both Ne vm 770 A (left) and Fe xn 195 A (right). Both light 
curves are obtained by averaging over 9" along the slit. Panels as in Fig.|4] 



